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Semiclassical  Theory  of  Inelaatic  Scattering  of  a  Particle 
in  the  Near- Adiabatic  Limit 


Lit-Deh  Chang  and  Walter  Kohn 

Department  of  Physics 
University  of  California 
Santa  Barbara  CA  93106 


ABSTRACT:  A  semiclaasical  theory  for  inelaatic  scattering  of  a  particle  by  a 
two-state  system  in  the  near-adiabatic  limit  is  developed.  The  exponentially  small 
transition  amplitude  is  calculated.  The  theory  is  an  extension  of  Pokrovskii  and 
Khalatnikov's  theory  for  above-barrier  reflection.  The  analysis  involves  studies 
of  the  WKB  solutions  in  the  complex  coordinate  plane  along  certain  contours. 
The  region  of  validity  of  the  theory  is  established.  Our  result  has  the  same  form 
as  the  Landau- Zener-S tuckelberg  formula;  however,  our  theory  is  applicable  to 
more  general  systems.  Numerical  comparisons  with  exact  solutions  are  presented; 

the  differences  between  our  results  and  exact  solutions  become  negligible  as  the 

i 

adiabatic  limit  is  approached. 
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L  INTRODUCTION 


The  inelastic  scattering  problem  of  a  particle  by  dynamical  systems  with  internal  de¬ 
grees  of  freedom  is  an  important  problem  in  many  fields  of  physics  and  chemistry.  In 
various  contexts  the  internal  degrees  of  freedom  could  represent  phonons  in  a  solid,  vi¬ 
brational  states  of  a  molecule,  atomic  energy  levels,  ionisation  states,  etc..  Consider  the 
scattering  of  a  particle  by  a  harmonic  oscillator  for  example.  In  the  adiabatic  limit  in 
which  a  heavy  particle  is  incident  with  low  velocity,  the  oscillator  wave  function  follows 
the  interaction  potential  due  to  the  particle  adiabatically  and  returns  to  its  initial  state  as 
the  particle  leaves  the  interaction  region;  even  if  inelastic  scattering  is  energetically  possi¬ 
ble,  the  inelastic  transition  amplitude  approaches  zero  in  the  adiabatic  limit.  For  a  small 
but  finite  incident  velocity  the  transition  amplitude  is  exponentially  small  as  a  function  of 
a  parameter  which  characterizes  the  deviation  from  the  adiabatic  limit.  Near  the  adiabatic 
limit  a  WKB  expansion  in  the  appropriate  small  parameter  is  natural,  such  an  expansion 
yields  a  vanishing  transition  amplitude  in  any  finite  order.  Many  approaches  1  have  been 
developed  in  the  past  to  deal  with  the  problem  with  various  degrees  of  success.  However, 
questions  of  the  validity  and  generality  of  those  approaches  have  remained  unresolved.  In 
this  work  we  have  developed  a  rigorous  treatment  of  this  problem  to  obtain  exponentially 
small  transition  amplitudes  for  general  interaction  potentials  between  a  particle  and  a 
system  with  two  internal  states. 

Studcelberg  2  was  the  first  to  applied  a  Zwaan  type  of  analysis  of  the  WKB  solutions 
to  the  problem  of  the  inelastic  transition.  His  theory  was  based  on  the  so-called  Stokes 
constant  method.  Since  the  Stokes  constant  method  has  not  been  justified  for  general 
problems,  the  validity  of  his  theory  remains  unclear  •*.  Our  theory  is  an  extension  of 
Pokrovskii  and  Khalatnikov’s  theory  *  for  above- barrier  reflection.  We  obtain  the  same 
result  as  the  Landau- Zener-Stuckelberg  formula  in  the  near-adiabatic  limit.  Our  theory 
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is  however  mors  general  and  ws  establish  its  region  of  validity.  In  Ssc.  II  ws  review  the 


work  of  Pokrovtkii  and  Khaiatnikov.  In  Ssc.  HI  ws  review  the  two-state  model  and  the 


adiabatic  basis.  Our  procedure  for  obtaining  the  transition  amplitude  is  presented  in  Sec. 


IV.  Numerical  results  are  described  in  Sec.  V.  Mathematical  proofs  of  the  propositions 


used  in  Sec.  IV  are  given  in  the  appendix. 


i  ' 


p  M  y  v  ;r. 


IL  ABOVE-BARRIER  REFLECTION 


Consider  a  one-dimensional  potential  scattering  problem  with  a  localized  repulsion 
potential  V (x).  IS  E  >  V{x)  everywhere  on  the  real  axis,  the  Schrodinger  equation 

+  a2 (B  -  K(*))*(*)  .  0  ,  (1) 

,  P) 

has  a  physical  solution  corresponding  to  incidence  from  the  left,  which  satisfies  the  following 
asymptotic  conditions: 


X  — *  — oo 


where 


q  =  lim  y/E  —  V  (x) 

x  *±oo 


R  is  the  reflection  coefficient  and  T  is  the  transmission  coefficient.  In  the  semidassical 
limit  a  -*  oo,  T  tends  to  one  and  R  is  exponentially  small  in  a.  On  the  real  axis  the  WKB 
solutions  of  Eq.(l)  are 

9 

±ia  f  q(x)dx 

7  *  (6) 


where 


q{x)  =  y/E-V(x) 


and  z  is  an  arbitrary  but  fixed  lower  bound  which  may  be  real  or  complex.  These  solutions 
have  an  error  of  order  a""1.  In  the  following  we  shall  use  the  symbol  &>  to  denote  an 
approximate  equality  which  becomes  exact  in  the  limit  a  — *  oo.  According  to  the  boundary 
condition  (3)  ,  the  physical  solution  can  be  expressed  in  terms  of  the  WKB  solution: 

¥(x)  » A+y9+{x,z)  ,x-+oo,  (8) 


4 


where  A+  is  a  constant.  As  we  trace  the  solution  along  the  real  axis  beginning  at  +oo, 
Eq.(8)  remains  valid  uniformly  all  the  way  to  —  oo.  But  does  not  include  the  exponen¬ 
tially  -<rwa.ll  reflected  wave.  In  fact  R  cannot  be  obtained  even  if  we  expand  the  solution 
to  higher  orders  in  powers  of  a-1  .  Pokrovskii  and  Khalatnikov  4  developed  a  method  for 
obtaining  the  exponentially  small  R  to  leading  order.  We  now  review  their  work  briefly. 

Throughout  thi*  work  we  shall  use  z  to  denote  a  complex  coordinate  and  use  x  to 
denote  a  real  coordinate.  We  assume  that  all  potentials  can  be  analytically  continued  into 
those  regions  of  the  z  plane  which  are  relevant  to  our  discussions.  In  order  to  obtain  R  we 
must  construct  a  solution  of  the  Schrodinger  equation  along  a  path  in  the  z  plane  where, 
even  as  a  — *  oo,  the  ratio  of  reflected  to  incident  wave  functions  (  see  Eq.(3)  )  remains 
finite  and  does  not  tend  exponentially  to  zero.  This  path  passes  through  a  point  Zc,  defined 
by  q(zc)  =  0  5,  where  the  WKB  solution  breaks  down.  However,  similar  to  the  familar  case 
of  a  classical  turning  point  on  the  real  axis,  near  Zc  an  exact  solution  of  the  Schrodinger 
equation  can  be  obtained  and  joined  to  WKB  solutions.  If  there  are  more  than  one  Zc,  the 
one  closest  to  the  real  axis  must  be  chosen. 


Consider  now  the  transition  point  z«.  The  appropriate  lines  passing  through  zc  are 
given  by  the  condition. 


Im 


z 


=  0 


(9) 


They  are  the  so-called  anti-Stokes  lines,  on  which  the  functions  tf+(z,  zc)  and  ¥~(z,zc) 
have  the  same  magnitude,  regardless  of  the  value  of  a  (  see  Eq.(6)  ).  There  are  altogether 
three  anti-Stokes  lines,  two  of  which  are  shown  schematically  in  Fig.  1.  L\  and  £>2  are 
useful  to  us  because  they  form  a  contour  connecting  +00  and  —00.  Asymptotically  they 
run  parallel  to  the  x-axis  at  a  distance  which  we  denote  by  yi .  The  procedure  of  Ref.  (4) 
for  obtaining  R  is  as  follows. 


Starting  with  the  right-going  wave  (4)  at  +00  on  L\ ,  we  find  the  wave  function  in  the 


5 


asymptotic  region, 


*M  -  ,  (io) 

where  z  =  x  +  tyj.  Except  in  a  small  neighborhood  of  radius  0(a-2/3)  near  2c,  'I'(z)  can 
be  represented  on  L\  by 


¥(2)  a  A+¥+(z,Ze) 


where 


and 


A+  =  re"ia,,+ 


±00 

■  ±{  J  (*(■*')  -  f)dz'  -  ?*c} 


(11) 


(12) 


(13) 


The  solution  is  then  matched  to  the  exact  solutions  of  the  Schrodinger  equation  in  the 
vicinity  of  Ze,  which  are  the  familiar  Airy  functions.  The  Airy  solutions  are  valid  in  a 
neighborhood  He  of  radius  dc  around  2^  (  see  Fig.  1  ),  where  dc  is  the  range  in  which 
the  term  {E  —V  (2))  can  be  approximated  by  the  linear  term  of  its  Taylor  expansion  at 
^ci  dc  is  evidently  independent  of  a.  The  matching  region  lies  between  |z  —  zc|  =  dc  and 
|z  —  ze|  ~  a-2/3,  a  must  be  large  enough  so  that  the  matching  region  exists.  Similarly, 
except  for  a  small  neighborhood  of  Ze,  the  physical  solution  on  Li  can  be  represented  by 
a  linear  combination  of  the  two  WKB  solutions,  Eq.(6), 


»(*)  »  {B+¥+(z,Zc)  +S-¥-(z,Zc)} 


(14) 


where  f?+  and  B~  are  determined  by  a  similar  matching  procedure  to  the  Airy  solution 
in  the  vicinity  of  2c.  B+  and  B~  have  the  same  magnitude.  On  Li,  as  z  — *■  -00,  ^(z), 
Eq.(14)  takes  the  form 


•(*) 


a 


-L{B+e~,a,1~  eiaqi  +  B~  eiav~  e~iaqz 
V? 


} 


(15) 


When  the  appropriate  values  of  B+  and  B  are  used  and  comparison  with  Eqs.  (3) 
and  (4)  is  made  one  finds 


T « eia(1’*+Tf~)  ,  (16) 

R « -ie2iav~  .  (17) 


Notice  that  |T|  as  1,  while  R  is  exponentially  small. 

Let  us  note  in  passing  that  R  cannot  be  obtained  by  using  a  contour  consisted  of  L* 
and  £>2,  corresponding  to  L\  and  L?,  but  passing  through  z*,  the  complex  conjugate  of 
Zc-  The  reason  is  that  whereas  on  L\  and  L% ,  the  asymptotic  magnitudes  of  the  incident, 
transmitted  and  reflected  waves  are  all  equal,  on  L*  and  ££,  the  magnitude  of  the  reflected 
wave  is  exponentially  smaller  than  those  of  the  incident  and  transmitted  waves  by  a  factor 
of  e~4a^n .  Such  an  exponentially  small  term  is  “missed"  by  our  procedure. 


m.  ADIABATIC  BASIS 


We  consider  the  one-dimensional  scattering  of  a  particle  of  mass  M  by  a  two-state 
system  with  energy  splitting  hu.  We  denote  the  internal  states  by  |0),  the  lower  energy 
state,  and  |1),  the  higher  energy  one.  They  form  the  so-called  diabatic  basis  for  a  two 
dimensional  vector  space.  The  Hamiltonian  of  the  total  system  can  be  written  as 


H  —  3P  +  3t  +  3i„t  * 

(18) 

*2  d2 

3p~  2Mdxz  + 

(19) 

hut 

Ha - - 

(20) 

Sint  =»  Vli*)*:  +  V2(x)ax 

(21) 

where  3P  is  the  Hamiltonian  of  the  particle,  B » is  the  Hamiltonian  of  the  two-state  system, 
H{ni  is  the  interaction  Hamiltonian  and  <7X,  <rr  are  Pauli  matrices.  The  two-component 
Schrodinger  equation  in  this  basis  is 

Ui(*)j  \*i(*)i 

We  assume  that  Vq(x),  Vi(x)  and  V^x)  all  tend  to  zero  as  x  -+  oo  and  that  the  particle 
comes  in  from  the  right.  For  a  fixed  energy,  the  velocity  of  the  incoming  particle  decreases 
as  the  mass  is  increased.  In  the  limit  of  large  a  the  particle  acts  on  the  two-state  system 
as  a  slowly  varying  external  field  and  the  two-state  system  follows  nearly  adiabatically;  in 
the  infinite  a  limit  there  is  no  inelastic  transition.  In  the  extreme  adiabatic  limit  (  a  =  oo 

)  the  Hamiltonian  is  trivially  diagonalized  at  each  z  by  the  transformation1,2 

Ceos$(z)  aind(z)  \ 

(23) 

— sintf(z)  cosO(z)  J 

where 

,<m  •  (24) 

with 


(25) 


vi(»)  *  vi(»)  -  y 

(25)  does  not  determine  9(z)  uniquely;  we  make  the  choice  of  0(z)  to  tend  to  0  as  x  — ►  oo 
so  that 


— w-4.  i— -jaa — 

A’W  +  v22W 


«■«<(*)  =  — t  (  1  +  ,  I  (27) 

slvt{.)+V*[m)) 

We  shall  assume  that  there  is  a  pair  of  points  Zc  and  z*  where  V\2  +  V 2 2  vanishes.  If 
there  axe  more  than  one  pair  of  we  chooee  the  one  which  is  closest  to  the  real  axis. 


y/V{*  +  Viz  is  defined  to  be  positive  along  the  real  axis  with  one  cut  going  from  z<-  to  10c 
and  the  other  from  z*  to  —too.  This  choice  of  branch  cuts  insures  that  the  mixing  angle 
9 (z)  is  a  continuous  function  of  x  on  the  real  axis.  For  z  real,  U  is  unitary.  The  new  basis 
vectors 


v  ii»*>  /  vn>y 

constitute  the  so-called  adiabatic  basis.  In  this  basis  the  Hamiltonian  takes  the  'om. 


Ha  =  UHU~ 


+  Ticrt*  )2  +  uo(*) 


*>*  all  1  ‘  a'  1 
TV9  ~  77 

jL  *  - 
-TWT^  ~  !T7  9  * 


where  Uq,U\  are  the  adiabatic  potentials. 


0o(*)  =  v<)(z)  -  v  V’:(*'r  - 


Ux(z)  =  Vn(z)  -  y,  V: 


•o(*) 

*i(*) 


The  Schrodinger  equation  becomes 


where 


It  is  easy  to  see  that 


E* 


f*o(*)\ 


=  E 


•oM  \  /  *q(z)  \ 

Ui  (z)J 


(33) 


(34) 


*0(z)\O)  +  »i(j)|1>  S  «o(*)|0,*>  +  *i(a)|lf  a)  (35) 

If  the  off-diagonal  terms  of  3a  are  neglected,  and  #i  are  uncoupled.  The  WKB 
solutions  (  a  oo  )  for  these  uncoupled  equations  have  the  usual  form, 

I 

±ia  f  q0(z)dz 

*(*,jM)"7SPr  7  ’  (36) 

and 

M 

±ia  f  qi(z)dz 

7  •  <"> 

where 

«(*)  *  y/E-Ui(z)  ,  (38) 

and  is  an  arbitrary  but  fixed  lower  bound.  The  branch  of  «(*)  is  defined  so  that 
9i(x)  —  y/\E  ~  VJ(*)|  for  x  >  Xj  and  ?,-(x)  =  —  Vi(x)|  for  x  <  zt-,  where  zt-  is  the 

i  1 

classical  turning  point  defined  by 

E  -  U{(xi )  =0  ,  i  =  0, 1  (39) 

If  we  treat  these  solutions  at  the  classical  turning  points  z,-,  we  will  obtain  a  reflection 
coefficient  of  magnitude  1  and  a  vanishing  inelastic  transition  amplitude.  This  remains 
true  in  all  higher  order  WKB  approximations.  We  need  to  go  into  the  complex  z  plane 
and  follow  the  solutions  along  appropriate  contours.  We  need  to  consider  the  singular 


points  :  xq,  x\  and  (  jm  Fig.  2).  Not*  that  at  *c  the  function  0(z)  diverges  and 
hence  the  off-diagonal  terms  of  Ha  also  diverge,  which  signals  the  failure  of  the  adiabatic 
solutions.  In  Fig.  2  we  show  the  anti-Stokes  lines  defined  as  follows, 


z 


LuL2  :  Im y  (flo(z)  -  ?i(z))  dx  *  0 

(40) 

X 

Lz  :  Im  I  7i(z)dx  =  0  , 

J 

(41) 

*1 

X 

L4  :  Im  f  qo(z)d*  =  0 

(42) 

* 0 


Li  and  L\  are  just  the  anti-Stokes  lines  for  ordinary  potential  scattering  problems  with 
classical  turning  points  at  x\  and  xo.  L\  and  L 2  are  the  lines  on  which  (z,  Zc)  ($q  (z,  Zc)  ) 
have  the  same  magnitude  as  $^"(z,  Zc)  ($~(z,  Ze)).  The  topology  of  Fig.  2  depends  on 
the  potentials  and  the  energy.  Since  U\  >  Uq,  x\  >  xq.  Both  xo  and  x\  shift  to  the  right 
as  the  energy  is  decreased,  whereas  *e  does  not  depend  on  the  energy.  In  most  previous 
discussions  of  Stuckelberg's  theory,  only  the  cases  where  Re(zc)  >  xq,  x\  are  discussed. 
We  will  see  that  this  restriction  is  not  necessary.  The  general  features  of  L\  and  L?  are: 
1.  L\  extends  to  infinity  and  becomes  parallel  to  the  real  axis  at  a  distance  yi  for  x  — *■  00; 
and  2.  Ly  intercepts  the  real  axis  at  a  point  2  which  can  easily  be  shown  to  be  always  less 
than  xi.  Both  y\  and  2  depend  on  the  energy  for  given  potentials. 


IV  INELASTIC  TRANSITION  AMPLITUDE 


Consider  the  physical  solution  of  the  Schrodinger  equation  (23)  along  the  real  axis  with 
the  following  asymptotic  conditions: 


»0(*)  -  e,a«*) 

-*  00 

-.0 

,x-»-oo 

(43) 

At 

,  X  — *  oo 

-*0 

,x  ->  -oo 

(44) 

where  $  =*  ?»(*  ~ *  oo).  In  the  adiabatic  limit  =*  0  and  |Aq"|  is  1.  In  the  near-adiabatic 
limit.  A*  cannot  be  determined  from  the  usual  treatment  of  the  WKB  solutions  to  all 
orders  whereas  Ajj*  can  be  determined  to  leading  order,  since  its  magnitude  is  —  1.  We 
have  seen  in  Sec.  m  that  at  the  complex  transition  point  Zc  the  adiabatic  solutions  break 
down.  We  shall  obtain  the  transition  from  one  state  to  the  other  by  a  careful  treatment  of 
the  Schrodinger  equation  near  s«:  we  need  to  obtain  the  exact  solutions  near  Zc  in  order 
to  connect  the  WKB  solutions  on  L\  and  Ij.  One  might  hope  to  obtain  the  exponentially 
small  amplitude  following  a  procedure  similar  to  that  described  in  Sec.  II;  however,  the 
geometry  of  L\  and  Ly  are  very  different  for  the  two  cases.  In  the  case  of  above-barrier 
reflection  Ly  extends  to  —  oo  where  WKB  solutions  become  the  exact  solutions  and  yield 
the  exponentially  small  reflection  amplitude.  In  the  present  case  goes  down  to  the  real 
axis  where  the  potentials  do  not  vanish  and  hence  the  exact  solutions  are  not  known.  A 
different  strategy  is  needed  in  making  use  of  the  solutions  on  L\  and  L^. 

It  is  helpful  to  think  of  the  physical  solution  as  a  linear  combination  of  two-component 
Joet  functions  (JF)  8.  A  JF  is  a  (  non-physical  )  solution  of  the  Schrodinger  equation 
which  asymptotically  has  only  one  incoming  or  one  outgoing  wave  in  one  channel.  The 
exact  physical  solution  corresponding  to  the  asymptotic  conditions  (42)  and  (43)  can  be 


written  as. 


*oM  -  +  Af*o*>(*)  +  ^T*oI+>(x>  •  («) 

*l(*)  =  +  A+’»<10+)(x)  +  A+*<1+,(x)  ,  (46) 

where  (x)  i a  the  component  in  channel  *  of  the  JF  with  unit  incoming  (— )  or  outgoing 

(+)  amplitude  in  channel  j.  Our  objective  is  to  determine  the  coefficient  A+  which,  by 
Eq.(44),  is  the  transition  amplitude.  The  procedure  is  as  follows. 

1.  We  first  obtain  A^,  whose  magnitude  is  —  1,  by  ordinary  WKB  techniques,  applied 
entirely  within  channel  0.  (  Ordinary  reflection  by  a  barrier.  ) 

2.  Next  we  note  that,  as  x  — *  —  oo  each  of  the  JF  components  in  Eqs.(45)  and  (46)  grows 
exponentially  like  4*,  Eq.(37),  which  is  the  most  rapidly  .growing  exponential.  We  shall 
choose  A*  so  that  the  total  coefficient  of  4*  vanishes. 

3.  To  accomplish  step  2,  for  Eq.(46),  we  need  to  know  the  coefficient  of  4+,  for  z  — ►  — oo, 
of  each  of  the  JF  components  4^°“^(z),  4^0+^(z)  and  4^1+^(x).  The  first  two  coefficients 
are  exponentially  small  and  need  to  be  carefully  calculated  on  appropriate  contours  (  see 
below  ).  The  last  coefficient  is  of  order  1  and  can  be  easily  obtained  by  ordinary  WKB 
techniques,  applied  entirely  to  channel  1. 

4.  Since  this  procedure  determines  A+  uniquely,  it  is  not  surprising  that,  in  fact,  it  also 
makes  the  leading  exponentially  growing  term  of  Eq.(45)  equal  to  zero. 

We  shall  need  a  number  of  propositions  which  are  proved  in  the  appendix  and  which 
have  been  verified  in  our  numerical  calculations. 

Proposition  1.  Consider  the  Schrodinger  equation  (33)  on  the  real  axis.  Under  con¬ 
ditions  detailed  in  the  appendix,  the  physical  solution,  which  vanishes  at  —  oo  and  has  an 
incoming  wave  in  channel  0,  has  the  following  components, 

*•!*>  *  (H 


) 


2 


A*(— a3£g(x))  ,  everywhere 


(47) 


»l(g) 


aO  ,x<*o. 


where 


*o(*) 

al^9\(x)**0  ,x>xq. 


X 

\  (^i(*))3/2  =  j  Vi** 


(48) 


(49) 


x* 


and  qi  is  given  by  Eq.(38).  The  branch  of  &  in  Eq.(49)  is  defined  so  that  £,-(z)  = 
[—U,(xi)\l^{z  —  xt)  in  the  vicinity  of  xt-.  Proposition  1  establishes  the  validity  of  the 
ordinary  WKB  solution. 

Proposition  2  gives  the  exact  solutions  in  the  vicinity  of  Zc- 

Proposition  2.  Consider  a  circle  f2e  of  radius  dc  around  z^.  We  denote  the  intersection 
of  flc  with  L\  and  by  *i  and  i«j  respectively  (see  Fig.2).  Assume  that  V{Vi  +  V^V 2  does 
not  vanish  at  Zc.  Then  for  sufficiently  small  dc,  the  general  solutions  of  the  Schrodinger 
equation  (23)  in  fle  are  given  by 

¥o(*)  »  (a+Aif— Ot2/3C2^?)  +  &+i?*'(-a2/3C2/3f)) 

+  (a_At(-a2/3C2/3f)  +  4_5*(-a2/3C2/3f))  ,  (50) 

to  the  leading  order  in  f,  where 

9c  =  -  vo{*c)  , 


C  = 


2?c 


(51) 

(52) 

(53) 


and  o±  and  6±  are  arbitrary  constants. 

Proposition  3  and  4  give  the  WKB  solutions  on  L\  and  Li  which  are  needed  for  finding 

*rw. 


Proposition  3.  Consider  the  anti-Stokes  line  L\.  Let  z\  denote  a  point  on  L[  at  a 
fixed  distance  from  Zc  and  let  L ^  denote  L\ ,  excluding  the  segment  from  zc  to  z\.  Under 


conditions  detailed  in  the  appendix,  there  is  a  solution  of  the  following  form, 


(54) 

(55) 


Proposition  4.  Consider  the  anti-Stokes  line  L?.  Let  z?  denote  a  point  on  L<i  at  a 
fixed  distance  from  Zc  and  let  LJ  denote  Lq,  excluding  the  segment  from  Zc  to  z%.  Under 
conditions  detailed  in  the  appendix,  there  is  a  solution  of  the  following  form. 


*0(*)  *  ,  (56) 

*l(»)  » *?(*,*)  .  (57) 

for  z  on  LJ  • 

On  Ly,  4g*  (z,  Zc)  and  (z,  z^)  have  the  same  magnitude.  Therefore  their  coefficients 
can  be  determined  accurately.  To  the  left  of  the  point  2,  (x,  Zc)  becomes  exponentially 

larger  than  $q"(x,Zc),  since  the  real  part  of  the  exponent  of  the  former  is  larger  than  that 
of  the  latter  for  x  <  2. 


We  now  proceed  to  find  the  leading  exponentially  growing  parts  (  for  z  — ♦  —  oo  )  of 
the  JF  components.  Consider  first  the  JF  component  ^g0+^(z)  in  Eq.(44).  As  i  — ►  +oo, 
^g0+^(z)  on  L\  is  an  outgoing  wave  in  channel  0, 

*<,0+)M  -  ,*-+00  ,  (58) 

where  z  =  x  +•  tyi .  According  to  proposition  3,  the  continuation  of  this  function  on  L\  is 
given  by 

*o°+)(z)  **  (*»*)«"*(*)  •  (59) 

where  $Q*(z,Ze)  is  the  WKB  approximation,  Eq.(37),  of  the  component  of  the  exact  solu¬ 
tion  along  the  adiabatic  basis  vector  |0,z); 


oo 

1c  =  J  {%-  (fo)dz  -  %Zc 


and  cas0(z)  is  given  by  Eq.(27).  The  approximation  is  valid  on  L\,  from  +oo  up  to  a  point 
*1,  where  |«i  -  *e|  =  0(a~2/ 3). 


We  now  match  this  WKB  solution  (59)  to  the  appropriate  combination  of  Airy 
functions,  Eq.(50),  at  a  point  at  which  both  the  WKB  solutions  and  the  Airy  solutions 
become  exact  in  the  limit  a  — *  oo.  By  taking,  for  example  |mi  —  Sc|  —  a-1/3,  we  see  that 
lml  ~  *c\J\*l  “  *c\  ~  a1/3,  validating  the  WKB  solutions,  and  |n»i  —  Sc|/|2i  —  *e\  ~~  a-1/3, 
validating  the  Airy  solutions.  At  mj  we  can  use  the  asymptotic  form  of  the  Airy  functions, 
since  the  argument,  |a2/3C2/3(mi  —  s«)|,  ia  ©(a1/3): 


2x  .  ,  2ir 

-  y  <  ar?(w)  <  y 


Accordingly,  inside  ne,  the  JF  component  is  given  by 


{A»(-o2/3C2/3f)  -  »S*(-cr2/3C2/3f)} 


where  C  is  given  by  Eq.(53)  and 


(65)  ; 


In  the  same  way  we  match  the  Airy-type  solution  (64)  inside  fJc,  to  the  appropriate  WKB 
solution  on  L?,  where  the  asymptotic  forms  of  Airy  functions  can  be  found  using  the 


following  identities: 


*•(,-*»/»„)  =  A.(w)  +  j«(»)} 


•Nfii 


♦y**, 


« 

'ft 


This  yields  ,  on  Z^, 

♦o°+)M  *  <>0+(,, *)<=<»«(*)  -  ^<l^)*i'(i,  Zc)nnl(z)  ,  (68) 

- c>  .  (69) 

where  t  =  ±1,  end  for  any  specific  potentials  is  given  by 

e  =  ^lim  ieatt(x)  (70) 

To  be  definite  we  choose  €  =*  +1  for  the  following  calculations.  From  Proposition  4,  Eq.(68) 
is  valid  on  the  whole  down  to  the  real  axis.  (  See  Fig.  2.  )  On  the  real  axis,  to  the 
left  of  the  point  2,  the  JF  component  is  dominated  by  (z,  Zc). 


Therefore,  on  the  real  axis,  we  have 

<+)M  - 


_^(°+)(x)J|n^(x) 


,z  —  +oo, 


,z  <  2 


where 


,  (73) 

A<$  =  e-<rre»an-ta^o  ^  (74) 

00 

m  =  Jfa(*)  -  h)dx  -  fcx,  ,  (75) 

and 

7  +  *7  *  J  -  J  q\{z)dx  (76) 

*0  *1 

It  is  easy  to  see  that  -j  is  positive,  hence,  is  exponentially  small.  The  0-component 
of  the  JF,  tfj,0  ^  (z),  needed  in  Eq.(44)  can  be  obtained  simply  by  taking  the  complex 
conjugate  of  ¥q0+^(x): 


<■’(*)  - 


«*  —  $(O-^(x)a*n0(x) 


,Z  —  +00, 


,z  <  2  , 


‘iS*!* 


i*. 


wi  w  *f(x’xi) 
=  "-<**“***+*** 


(79) 

(80) 


Finally  we  require,  for  Eq.(44),  #g1+\  the  0-component  of  the  JF  corresponding  to 
channel  1  outgoing.  In  the  adiabatic  basis  this  JF  is,  for  our  purpose,  obtainable  every¬ 
where  on  the  real  axis  by  ordinary  WKB  turning  point  methods: 

*(1+)r  '  1 


‘  [  ,  *  -let til  X 

ri  '  y/E 


,x  — *■  +00 

,X  <  Xl 


where 


and 


j(1+)  m  e~iarh 
Al,*i  e 


*^M»o 
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(1+) 


(*) 


,  everywhere  . 


(81) 

(82) 

(83) 

(84) 


The  reason  is  that  for  this  JF  a  single  component,  #[1+^(x),  is  dominant  everywhere  on 
the  x-axis;  this  is  aoi  the  case  for  the  JF’s  (0+)  and  (0—),  for  which  $q0±^  (x)  is  dominant 
for  x  >  X,  while  i^0±^(x)  is  dominant  for  x  <  i.  The  required  component  ^g1+^(x)  in  the 
diabatic  basis  is  given  by 

'®q1+^(z)  »  — #j1+^(x)«w*0(x)  ,  everywhere  (85) 

We  are  now  ready  to  determine  the  required  coefficient  A~f  from  the  condition  that  the 
coefficient  of  the  dominant  exponential,  for  x  — ►  — oo,  in  Eq.(44)  vanishes.  Substituting 
Eqs.(72),  (78)  and  (85)  into  Eq.(44)  gives, 


4.  a+  )  1  —  g 

A\,xx  +  *0  *l,xi  ^  Al  Al,xi  ~  u 


(86) 


For  a  wave  incident  in  channel  0  with  unit  amplitude. 


Substituting  in  Eq.(86)  for  Aq\  A^”\  A^^  and  A^j\  we  obtain,  on  solving  for  A+, 

Af  m  -2sm(o»7)e”aT'e,'a<r'0+^)  .  (88) 

Incidently,  we  can  verify  that  the  coefficient  of  (x,  xi)  in  Eq.(45),  for  x  <  £,  also  vanishes 
as  it  should  if  A*  is  given  by  Eq.(88). 

following  the  standard  definition  of  the  S  matrix, 

*0<*)  -  -L  («-“«*-  So0‘”*'1)  ,*-+«  ,  (89) 

*l(x)  .X  — •  4-oo  ,  (90) 

we  find  that 

Sqo  m  ie2***  ,  (91) 

and 

Sqi  *  2 «»(an)e“an'e,a(T'0+^)  (92) 

( If  e  is  chosen  to  be  —1,  the  right  hand  side  of  Eq.(92)  would  change  sign.  ) 

Eq.(92)  has  the  same  form  aa  the  Laadan-Zener-Studcelberg  formula  in  the  near- 
adiabatic  limit.  However  whereas  the  earlier  derivation  was  limited  to  the  case  where 
Zc  is  near  the  real  axis  (  almost  crossing  adiabsts  on  the  real  axis  )  and  well  to  the  right  of 
the  classical  turning  points  xq  and  x\%  our  derivation  is  not  limited  by  these  restrictions. 

We  shall  now  indicate  the  conditions  under  which  our  near-adiabatic  approximation  is 
valid.  It  requires  the1  existence  of  admissible  matching  points  m^,  on  L\  and  m2,  on  L<i 
(see  Fig.2).  It  is  easy  to  see  that  if  rn\  exists  so  does  m2.  The  existence  of  m\  requires 
that 

|»1  -  *e\  <  <*e  •  (93) 

Here  z\  is  given  by  the  following  condition 


which  is  the  analog  of  of  the  condition  signalling  the  breakdown  of  the  standard  WKB 
approximation  7, 


1  d  1 


=  1 


a  dz  q(z) 

and  dc  is  the  distance  from  Zc  over  which  the  following  expansion  is  valid 


(95) 


»(»)  3  ^-V0W  +  V/Vi(i)2  +  V2(z)2 
=  He  +  c{j  -  A,)1/2  > 


(«•) 


where  and  C  are  defined  in  Eqs.(52)  and  (53).  Combining  (93)  and  (94)  gives 

o*4/2|C|  >  1  (97) 

Since  a  =■  ■s /M/K  it  is  evident  that  for  any  given  potentials  and  given  incident  energy,  the 
condition  (97)  will  always  be  satisfied  for  sufficiently  large  M. 


The  left  hand  side  of  Eq.(97)  is  the  dimensionless  large  parameter  which  characterizes 
the  deviation  from  the  adiabatic  limit.  (  Let  us  remark  that  condition  (97)  is  not  necessarily 
equivalent  to  the  intuitively  suggestive  criterion  that  the  exponent  cry,  in  Eq.(92),  be  :»  1.) 


V  COMPARISON  WITH  NUMERICAL  RESULTS 


We  present  comparisons  between  our  results  and  exact  numerical  solutions.  We  choose 
exponential  potentials, 


Vb(z)  =  uoe~* 

t 

(98) 

vl(*)  -  “ 

flu 

'  T 

(99) 

V2(z)  -  u2«-* 

f 

(100) 

where  the  are  the  strengths  of  V)(a).  Varying  u,-’s  we  can  study  both  curve  crossing  and 
curve  noncroesing  cases.  We  shall  set  tiu  =  1. 

1.  Curve  non-crossing.  In  thia  case  we  choose  uq  =  1.02,  uj  =  .02,  U2  =  .2.  Two 
energies  JE?  =  1, 10  are  studied  for  different  values  of  a.  The  two  cases  have  very  different 
topology  of  *e,  xq,  x\  and  anti-Stokes  lines  as  shown  in  Figs.  2(a)  and  2(b).  Although 
the  diabatic  potentials  Vq(x)  ±  Vi(x)  cross  at  x  —  -3.2,  this  point  is  far  to  the  left  of 
the  classical  turning  points.  In  the  region  where  the  transition  takes  place  the  adiabatic 
potentials  are  almost  parallel  to  each  other.  Du  Fig.  3  the  absolute  values  of  $q0~^(x) 
and  *(10+)(x)  for  E  »  10  are  plotted  for  different  values  of  a.  For  large  at,  channel  1  is 
seen  to  dominate  for  x  <  2.  Numerical  results  for  the  S  matrix  are  listed  in  Table  I.  The 
differences  between  semiclassical  and  exact  solutions  diminish  in  the  limit  of  large  a  for 
both  energies. 


2.  Curve  crossing.  In  this  case  we  choose  uq  =  1.2,  =  .5,  u2  =  .2.  The  diabatic 

potentials  cross  at  x  =  0.  Two  energies,  E  =  0.7, 1.5,  are  studied.  For  E  =  1.5,  the 
topology  of  2c, xo,  xi  and  anti-Stokes  lines  is  similar  to  Fig.  2(a)  and  for  E  =  0.7,  it  is 
similar  to  Fig.  (2b).  Results  of  the  S  matrix  are  also  included  in  Table  I.  The  differences 
between  semiclassical  and  exact  solutions  also  diminish  in  the  limit  of  large  a  for  both 
energies. 
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VI  CONCLUSION 


We  have  presented  a  rigorous  treatment  of  the  semiclassical  theory  for  obtaining  ex¬ 
ponentially  small  inelastic  transition  amplitudes.  The  region  of  validity  of  our  theory  is 
established.  Our  result  ,Eq.  (92),  has  the  same  form  as  the  Landau-  Zener-S  tucJceiberg  for¬ 
mula  in  the  near-adiabatic  limit.  However,  our  theory  is  more  general  and  is  applicable  to 
a  wide  range  of  systems. 
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APPENDIX 


In  the  following  proposition*  we  make  the  same  assumptions  about  the  potentials  as 
in  the  text.  The  proofs  of  propositions  1,  3  and  4  are  analogous  to  Jeffreys’  proof  8  for  the 
ordinary  potential  scattering  problem.  For  convenience  we  shall  use  N  to  denote  a  bounded 
constant  independent  of  a.  The  value  and  the  dimensionality  of  N  may  be  different  at 
different  places  and  are  irrei event  to  our  discussions. 


I.  Proposition  1.  Consider  the  Schrodinger  equation  (33)  on  the  real  axis.  We  assume 
that  for  the  given  energy  E  each  adiabatic  potential  U{,  i  =  0, 1,  has  only  one  classical 

turning  point,  and  that  U-  (xt)  <  0,  which  will  normally  be  the  case.  If  the  following 
functions  are  bounded, 

*M-/(lJtoMI  +  \MW)  /fgl) 

and  if  the  following  quantities  can  be  bounded  by  a1/6  IV, 

F*±  =  [  +  7\j_  f  y/fzXzQ _ 

JL  {  l<fa'  VV&<*')«i(*')(*(*')±» ,(!-))  J\ax' 

*» 

(*>j)  =  (0,1),  (1,0),  (103) 

;  I 

where  •  —  x,  ±  Na  2^3,  £,  is  given  by  Eq.(49),  /l;,  and  h  are  given  by  Eqs.(l09)-(112), 
then  the  physical  solution,  which  vanishes  at  -oo  and  has  an  incoming  wave  in  channel  0, 
has  the  following  form. 
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dx1 
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(101) 
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*o(*)  23  jj.  'j  A*(  —  a3fg(x))  everywhere  . 
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aI/6«1(x)»0  x  >  Xn , 


1104) 


MOM 


Proof.  Defining  a  new  set  of  variables 


x«(&(*))  -  {j£)  * ’ *  ~  °» 1 

(106) 

the  Schrodinger  equation  (33)  can  be  written  as 

^“?X0  +  <*2Soxo  *  /oo(£o)xo  +  /oi(£o)xi  + 

(107) 

^JXI  +  a2llXl  =  /ll(£l)xi  +  /io(£i)xo  -  M^l)^X0 

(108) 

where 

L/+\  2^ 

M£i)  =  (€o'ei')1/2  • 

MM  2(jj]J  4  «j)«  +  (®2  ’  ’  • 

(109) 

(110) 

/0l(&)  -  (-a,^,)3/;  (<"«l'  -  *(l")  , 

(111) 

Aottl)  *  (^/fl/)3/2  (  *  &  +  ) 

(112) 

In  Eqs.(l07)  and  (108)  x«  are  functions  of  &.  The  correspondence  between  fo  and  is 
one-to-one.  If  xi  i»  negligible,  Eq.(107)  involves  only  channel  0.  The  solution  of  xo  13 
given  by  Eq.(104).  To  show  that  xi  i»  indeed  negligible  we  transform  Eq.(108)  into  an 
integral  equation  8.  We  need  to  construct  the  Green’s  function.  The  general  solutions  of 
the  homogeneous  equation, 

^X(i0)  +  a2£lX(i0)  =  0  ,  (113) 

are  Airy  functions: 

*f>(«l)  =  ,  (114) 

where  ai  and  b\  are  arbitrary  constants.  The  required  Green’s  function  which  satisfies  the 
differential  equation, 

j2 

+  «2£iG,(£i,£i)  =  <S(£i  -  ?i)  ,  (115) 

<K1 
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has  the  following  form: 


Cl  (ft,  ft)  =  A.(-a2/3ft)iMft)  ,ft  <  ft, 

=  {-HA.(-<,V3ft)  +  fl.(-a2/3ft)}Pl(ft)  ,ft>ft,  (118) 
where  G\  — *■  0  as  £1  — *■  —  oo  and,  as  £1  — ►  +oo,  is  an  outgoing  wave  in  channel  1  for 

a* 

fixed  Solving  for  G\,  we  obtain, 

Gi(ft,  ft)  =■  -T5n^'(-«2/3ftK*^(-«2/3ft)  +  S.'(-o2/3ft)}  ,  ft  <  ft, 

A.(-a2/3ft){U»(-a2/3ft)  +  fl.(-a2/3ft)}  ,  ft  >  ft.  (117) 


a2/3 

T 

'^273 


The  integral  equation  for  the  required  solutions  is, 

+00 


Xl(£l)  *  J  ^?l<»l(£l.  ll)  j/llXl  +  /10XO  -  A^"Xo} 

—00 


(118) 


where  the  homogeneous  term  is  chosen  to  be  0.  (  We  have  omitted  some  arguments  of 
the  funcitans  to  shorten  the  notation.  )  The-  first  intention  of  Eq.(118)  is  obtained  by 
substituting  xi  with  0  and  xo  with  A»(— a2/3^): 

Xl(€l)  =  [  <<ftGl(ft,ft){/l0-*4-}Ai(-a2/3a)  (119) 

J  at  1 


—00 


The  right  hand  side  of  Eq.(119)  can  be  written  as  a  sum  of  three  terms, 

+00 

X^’tft)  -  .^jM-«2/3ft)  /  <<ftA-(-«2/3ftH/lO  -  /.-i}^,(-a2«&)  ,(120) 


-OO 


xf’tft)  =  -^73  J  <*ftiC(ft,ftK/io  -  *.4-}A.(-a2« 


£0) 


—00 

ex 


xf’tft)  =  ^73  /  <*ftJ^(ft, ft){/lO  -  *4-}Ai(-a2'3{0) 

+00  ^ 


(121) 


(122) 


where 
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-  B«(-a2/3e1)i4«(-a2/56) 

5)  -  M~ a2/3&)S»'(-a2/3S) 


(123) 

(124) 


Since  Ai(— a2/3^)  is  bounded  everywhere  and  A«(— g^/3£i)  is  exponentially  smaller  than 
Ai(— o2^3^))  for  x  <  xq,  we  have 


Xl  (x) 


S3  0  ,X  <  XQ 


A'(-« »/»&(*)) 

*0  ,x  >  xq 


(125) 

(126) 


Both  xi2)  end  xf  satisfy  equations  similar  to  (125)  and  (126).  We  shall  give  the  bound 


(2) 

for  the  second  term  of  x)  • 


xf(«l)  =  /  <£iCjv4-A.(-<»;V3&) 


(12T) 


the  bound  for  other  terms  in  x^  and  x^  can  be  obtained  similarly.  Consider  first  the 
region  x  <  xq.  In  order  to  utilise  the  asymptotic  form  of  A*'(-cr2/3£o)  we  disscuss  two 
cases  ,  x  <  and  x^  <  x  <  xq,  separately. 


For  x  <  Xq  ,  using  the  asymptotic  forms  of  Airy  functions,  we  have, 

.  or  /  \qildx"  f  -a(f  \qi\dx"+f  \q0\dx") 

{  /*'•«*<**  "  - 


where 


*w  -  -7^4 


V?i(x)^o(*) 


The  integral  in  Eq.(134)  can  be  written  as 


f  jotxffixQ  d  ---(/‘if. vrJfM 

]  -io(«o(i')+n(i'))iii'  ' 


(128) 


(129) 


Performing  an  integration  by  part  in  this  form,  we  obtain, 


am.  < ,  i  **(*) 

1  '  «7/6  (^  (*)  +  *!  (*))  1 

+  1  — -1  -  /  4*'(JL-J0£-)e  *  *  *  \  .(130) 

a7/6  ^i(x)  J  %  +  <n 

— OO 

Since  we  have  assumed  that  there  is  only  one  classical  turning  point  for  each  adiabaitc  po- 

*t  *0 

tential,  the  function  /  |$i| dx//  +  /  \%\da/'  is  monotonically  increasing  as  x  moves  from  xo  to 

X  X 

r° 

— oo.  Hence  the  exponential  term  in  the  integral  of  Eq.(130)  is  bounded  by  exp(— a  /  \qQ\dx") 

x 

From  this  and  the  assumption  about  F9+,  it  is  not  difficult  to  see  that 


lxi2)l  <  ~a*(-«2/3&(*)) 


t*  <  xn 


(131) 


For  Xq  <  x  <  xq  the  asymptotic  forms  of  Airy  functions  are  not  valid.  However,  both 
At  and  Bi  are  bounded,  we  can  estimate  the  magnitude  of  the  integral  easily.  In  this 
region  the  integrand  can  be  bounded  as  follows: 


<  £  . 


(132) 


-Ai(-a2/3f0)l  <  a2/3.'V 


(133) 


v  a 

Since  the  integration  region  is  only  of  0(a~ ,  the  integral  is  only  0(a-1).  Therefore 


we  conclude  that 


xf’W 

*o(*) 


ss  0  ,  x  <  XQ 


(134) 


For  x  >  xq,  Xj2^  can  be  analised  for  regions  xg  <  x  <  x1  ,  x1  <  x  <  x^  and  x  >  x* 
by  a  similar  treatment.  We  obtain  «  0. 

II.  Proposition  2.  Consider  a  circle  Qc  of  radius  dc  around  zc.  Assume  that  V[V\-{ -V^Vy 
does  not  vanish  at  Zc  then  for  sufficiently  small  dc  and  sufficiently  large  a,  the  general 


I 


solutions  of  the  Schrodinger  equation  (23)  in  Hc  is 


*0 (*)  »  (o+A»(a2/3C2/3f)  +  4+B*(a2/3C2/3f)) 

+  e-****  (a _^*(a2/3C2/3f)  +  *_S*(a2/3C2/3f))  ,  (135) 

to  the  leading  order  in  f,  where  f,  qc  and  C  are  given  by  Eqs.(51),  (52)  and  (53);  and  a± 
and  b±  are  arbitrary  constants. 

Proof.  Eliminating  ®i(s)  from  the  Schrddinger  equation  we  obtain  the  fourth  order 
differential  equation  for  ¥(](*), 

D*oW  *  {^  +  /jW^  +  /sW  Jr  +  /iWs  +  /ow}  *oW 

=  0  ,  (136) 

where 

fz 38  ~27^  »  (137) 

ft-  2(B-V0)+a2(-S+J&2)  ,  (138) 

k2  ^2 

h  =  -2a2(V0'  +  V,'  +  -  V0  -  n)|  ,  (139) 

/0  =  a \(E  -  V0)2V ?  -  V22) 

+  a2(-(v0"  +  V")  +  (g  -  r0  -  vt)  +  2^(V0'  +  y')|  .  (ho) 

Pulling  out  a  fast  oscillating  part  of  ¥o> 

*oM=-«i“,*<x(r)  ,  (i«) 

the  solutions  of  x(f)  include  two  slowly  varying  solutions  (  for  sufficiently  small  |f|)  and 
two  fast  oscillating  solutions;  we  will  find  the  former.  The  equation  for  x  reads 

I^  +  J2WS3 +I2(z)S2  +  »iw4  +  joM}xM  =o  -  (14J) 


where 


03(e)  *  —Aiaufe 

(143) 

92  (*)  *  2a2  {E  -  V0  -  3 q2) 

91  (*)  »  4o3wc(^  -  Vq  -  9?) 

(144) 

+  2a2[V0'  -  V'  -  (V0  -  Vj  + 

(145) 

*>(*)  *  »4[(fc4  -  -  V0)  +  (*-  V0)2]  -  (V2  +  V22)] 

+  -  VjJ  +  [E  -  V0  -  V,  -  4)^1 

(146) 

Linearizing  the  potentials  at  *c  and  changing  variable  from  *  to  f  we  obtain 

{£ 

+  °^3 +  («2^2  +  a3^2f)^J  +  a2Sl^  +  +  a4T0f)}  x(f) 

=  0 

1 

(147) 

where 

53--2g  . 

(148) 

*  -4q%  , 

(149) 

= -2ax  -  2oq  +  2^  +  4^^ 

®2  ®2 

(150) 

T2  =  4*0cao 

(151) 

50  =  *(fti<*2  -  «1*2) 

(152) 

To  =  2(0565  +  » 

(153) 

with 

of  =  v;'(zc)  , 

(154) 

*1,2  =  ^l,2(*c) 

(155) 

Making  a  scale  tranaformatin 


I  =  a2/3? 


1156) 


and  assuming  that  arc  of  the  »mm  order  as  |x|,  Eq.(153)  can  be  reduced  to  the 

leading  order  in  a: 

s?+cS,x+^^|530+r‘lr'SMl“0  (l57) 

For  tj  >  0( a~1/3)  the  last  term  is  negligible,  the  solutions  are  the  Airy  functions.  For 
H  <  0(a-1/3)  the  last  term  is  not  smaller  than  the  second  term,  however,  all  terms  are 
small;  the  correction  to  the  Airy  functions  vanishes  in  the  large  a  limit.  Hence  the  slowly 
varying  solutions  are 

X(f)  w  a+A»(a2/3C72/3f)  +  6-,.Ht(a2/3C2/3f)  ,  (158) 

where  &+  and  4+  are  arbitrary  constants.  This  shows  that  the  scaling,  Eq.(156),  is  indeed 
a  correct  one.  Eqs.(141)  and  (158)  give  one  set  of  solutions.  Pulling  out  another  fast 
oscillating  factor  e“ia^  (cf.  Eq.(141))  we  can  obtain  the  other  set  of  solutions  in  Eq.(135). 

HL  Proposition  3.  Consider  the  anti-Stokes  line  L\.  Let  z\  denote  a  point  at  a 

fixed  distance  from  %c  and  let  denote  L\,  excluding  the  segment  from  z<;  to  z\.  We 

assume  the  following  conditions  to  be  true:  1.  L\  does  not  pass  through  singularities  of 

z 

the  potentials.  2.  the  mappings  &(*)  =**  / <n(z)ds,i  »  0,1,  are  one  to  one  9;  3.  Im(£,) 
increases  monotonically  as  *  moves  to  the  right  on  L\  10.  If  the  following  functions  can  be 
bounded  by  N, 


00 

Fo(‘)  - 

J  (l/oo(*)l  +  l/oiWI)l«W^I  . 

(159) 

z 

OQ 

Fi(i)  - 

J  (l/ll(*)l  +  l/lo(*)l)  kl(*)da| 

(160) 

F?(t)  - 

(161) 

Z 


where  fij  and  h  are  given  by  Eqs.(187)-(170),  then,  for  z  on  L there  is  a  solution  of  the 
following  form, 

*o(»)  *0  ’  (162) 
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(163) 


ilif) 

*o(*) 


*  0 


where  ♦g'(i,  Zc)  is  defined  in  Eq.(36). 


Proof.  Defining  a  new  set  of  variables 


*(&(*) ) -(H)’ ♦■•<*)  .*'  =  0.1  • 

(164) 

the  Schrodinger  equation  (33)  can  be  written  as 

+  <*2xo  *  /oo(£o)xo  +  /oi(£o)xi  +  h(^o)^Xl 

(165) 

d2  d 

^-5X1  +  <*2Xl  *  /u(£i)xi  +  /lo(£l)xo  - 

(166) 

where 

Lt*  ^  20 

Mf,) "  (a' ft')1'5  " 

/«(«.)  4  (£)4  +  ({-)2  ■ 

(167) 

(168) 

/0l(&)  -  (&,{i,)3/2  {*  ft  *  ''ft  ) 

(169) 

/10  (ft)  =  (io,^,)V3  (-'"a'  +  *V)  . 

(170) 

The  correspondence  between  variables  £q  and  £1  is  one-to-one.  Similar  to  the  proof  of 
proposition  1,  we  transform  (165)  and  (166)  into  coupled  integral  equations  6: 


xo(£o)  =  Xo0)(^o)  +  j  ttoKatt o*fo)x 

j/oo(£o)xo(?o)  +  /oiUo)xi(ii)  +  ’ 

xiUi)»x(i0)Ui)-h  /  4iJT«(ft,?i)x 

{/ii(?i)xi(?i)  +  /io(£i)xo(£o)  -  MfO^-xidi) | . 


(171) 


(172) 
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where  Xg  '  end  X  [  ^  ere  solutions  of  the  homogeneous  equetions: 

Xo0)(&)  =  “Ot^0  +  4oe'“6  . 

X(i0)(fl)  =  «!«“*■  +  »ie-“<‘  , 


(173) 

(174) 


end  Ka  is  the  kernel, 

*.(&.  I)  =  - 1  -  «-"»« i-ZA  (ITS) 

Choosing  the  lower  bounds  of  the  integrals  to  be  ^°,  where  *(£f°)  =  oo+*yj,  and  choosing 

Xq0)  * 


x[0)»0 


(176) 

(177) 


the  integral  equations  become, 


Xo(6))  *  «*°*°  +  f  d£oKa  |  fooxo  +  /oiXl  +  h~xi  1  .  (178) 

£  1  ^0  J 

fi 

Xl (^l)  »  J  dtlKa  |/nXl  +  AoXO  -  ^T*0}  (179) 


The  first  iteration  of  Eqs.(178)  and  (179)  gives 

So  ^ 

Xotfo)  =  +  J  dtQKafweiat°  , 

*o°° 
fi 

Xl(*l)  -  J  d&Kc  (/10  -  h4r]  eiat° 

d*X 

Eq.(lSO)  includes  only  one  channel,  which  is  similar  to  the  ordinary  potential  scattering 
problem;  the  integral  is  negligible.  We  will  only  show  the  the  bound  for  the  second  term 
of  (181); 


(180) 


(181) 


si 

X^tfl)  5  -  f  dtxKah-i-joS* 
An  ®£l 


(182) 
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Performing  an  integration  by  part  (  which  is  similar  to  the  step  from  Eq.(l28)  to  Eq.(130) 
),  we  obtain 

X*!1’  -  (B,  +  B2  +  fl3  +  B„)  ,  (183) 

ot 

where 


p  =  Mf)gg(f) 

1  »(*)-«(*)  ’ 

T3  _ 

$o(*)  +  <ni*) 

BS  -  /  , 

y  Wflofso-aOiT 

oo 

P.  ,  f  dj(  d 

4  J  **[ <*«>(*)  + <nVY 


(184) 

(185) 

(186) 
(187) 


oo 

From  the  definition  of  ,  |e*a^°“^1)|  is  1.  Therefore  B\,  B%  and  £3  are  bounded  for  z  on 
.  From  condition  3  about  L\  we  see  that  the  real  part  of  the  exponent  of  £4  is  always 
less  than  0.  Thus  £4  is  also  bounded.  Therefore  we  have 


5&,0 

*oM 


(188) 


The  proof  of  proposition  4  is  similar  to  what  we  have  just  shown.  We  only  give  the 
statement. 


Proposition  4.  Consider  the  anti-Stokes  lines  L^.  Let  z<i  denote  a  point  at  a  fixed 

distance  from  Zc  and  let  denote  Ly,  excluding  the  segment  from  Zc  to  z 2.  We  assume 

the  following  conditions  to  be  true:  1.  Ly  does  not  pass  through  singularities  of  the 

z 

potentials.  2.  the  mappings  £,-(z)  =  f  qi(z)dz,  i  =  0,1,  are  one  to  one;  3.  Im(({)  increases 
monotonically  as  z  moves  down  to  2  on  11 .  If  the  following  functions  can  be  bounded 


by  N, 

Z 

Fq(z)  =  J  (I/00MI  +  |/oi(*)|)  Mz)dz\ 


(189) 


J 


(190) 


Z 

AM  =  J  (l/nMI  +  1/ioMI)  |giM<i*| 

*x 

p*±c-\  _  f  I  ^  f  ?iM)Mz/)  A 

(*W  k  Um)±«m); 


<M  ,*  =  o,i 


(191) 


where  /,,,  and  h  are  given  by  Eqs.(167)-(170),  then,  for  z  on  L^,  there  is  a  solution  of  the 
following  form. 


*oM  »  *o"(*,*e) 

where  $q(z,  Ze)  and  $i(z,  Zc)  are  defined  in  Eqs.(36)  and  (37)  respectively. 


(192) 

(193) 


£ 


r«i 

y> 


& 
f  ! 


S 
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FIGURE  CAPTIONS 


Fig.  1.  The  typical  topology  of  the  complex  transition  point  Zc  and  anti-Stokes  lines 
L[  and  are  shown  for  above  barrier  reflection  problem.  L\  extends  to  +00  and  becomes 
parallel  to  the  real  axis  at  a  distance  yi  asymptotically.  Li  behaves  similarly  as  it  extends 
to  —00. 

Fig.  2.  The  transition  points  zq,  *i,  2,  Zc  and  anti-Stokes  lines  L\,  Li,  £3,  £4  are 
shown  for  the  carve  noncrossing  case  in  Sec.  V:  2(a)  for  E  =  10  and  2(b)  for  E  =  1. 

Fig.  3.  The  absolute  values  of  9 jj0"^  (z)  and  9^^  (z)  for  the  case  of  Fig.  2(a)  are 
shown  for  several  values  of  a:  3(a)  for  a  =  1,  3(b)  for  a  =  5  and  3(c)  for  a  =  10.  |$q0+^(z)| 
is  represented  by  the  solid  line  and  |9^0+^(z)|  by  the  dash  line.  The  vertical  dot-dash  line 
marks  the  position  of  2.  Note  that  for  large  a,  $ (z)  becomes  greater  than  (z) 
as  x  moves  to  the  left  of  2. 

TABLE  CAPTION 

Table  L  The  numerical  results  of  exact  and  semiclassicai  solutions  for  the  S  matrix 
are  shown.  The  semiclassicai  results  are  calculated  from  Eqs.(9l)  and  (92).  Note  that  the 
differences  between  them  decrease  as  a  is  increased.  1  s 


Fig.  2C“) 


TABLE  I 


